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ABSTRACT 


The objective of this investigation is to determine the upper limits on the 
thermal yield of a l Kt source detonated at shallow depths. 
'I'hree problems 
were computed for the depths of burial of 0.5, 1.5 and 3 m. 
The source is 
an iron bubble with 30 an radius and 25 Kg of mass. 
The computation was 
done on the CEL code. 
The results show that the relative thermal yield is 
no more than 6%, and that it vanishes at OOB of about 4 m. 
.The computed 
blast strength is plotted against the range. 
The results are consistent 
with the measurements on "Johnnie Boy" and "Jangle-U", and in good agree- · 
ment with the RAD9 calculations and the HULL calculations. 


Introduction 


This report describes the results of a study on the thermal yield from 


the subsurface nuclear explosions. 
The objective of the investigation is to 


determine the upper limits on the thermal yield from a l Kt source detonated 


at shallow depths. 
The depths of burial (OOB) are 0.5, 1.5 and 3 m. 
The 
1,2 
computation was done on a Coupled-Eulerian Lagrangian (CEL) code. 


An explosion in the free atmosphere loses as much as 30 to 40% of its 
3,4 
energy by radiative processes such as x-ray emission and thermal pulse. 


Of these processes, more than. 95% of the radiative loss is due to thermal radi- 


ation. 
When the detonation takes place at shallow depths, the radiative process 


is inhibited by the optically· opaque vaporized soil and other high Z materials. 


The thermal yield rapidly decreases with increasing depth of burial. 


The peak over-pressure of the blast wave is plotted and compared with the 


other computations, the RAD9 calculation and the HULL calculation, and the measure- 


ments on "Johnnie Boy and Jangle-U". 
The CEL computation is consistent with the 


measurements. 
Within the computational error these three sets of ccalculations are 


mutually consistent. 
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Computational Technique 


The computation was done on CEL using the LLL equation of state and opacity 
5 
1,2 
tables. 
CEL is a two-dimensional -radiation - hydro code. 
It calculates 


the equilibrium temperature with a diffusion equation. 
CEL couples the Eulerian 


and Lagrangian grids, utilizing the advantages of both schemes. 
The Lagrangian 


formalism is iriherently more accurate than the Eulerian scheme for resolution of 
6 
shock variables. 
However, when there is a strong mixing of materials, the use 


of the Eulerian grid alleviates many computational difficulties. 
Therefore, 


CEL is most suitable where there is a strong interaction of gaseous materials 


with solids. 


Using an iron bubble of l Kt energy as a source, the explosions at 0.5, 1.5 


and 3 m were examined. 
Figure l shows the configuration of the problems. 
The 


atmosphere is represented by a Lagrangian grid because the accurate calculation 


of energy and other shock variables is required in this region. 
Since the ground 


motion is of little interest at this time, it is mapped in an Eulerian mesh. 


This alleviates computational difficulties at the expense of accuracy. 
The 


source is given by a Lagrangian grid so that an accurate account of the energy 


can be made. 
The sour.CP radius is 30 em, and the total mass is 25 Kg. 
The 


tracer line is an air-soil interface. 


These shallow explosions quickly communicate with the ground surface, and 


the sourre moves in an up.vard direction resulting in a "mushroom effect" (see 


Figure 2). 
'Ihe airblast appears as if it originated from a virtual point of 


explosion located at about 20 m above the ground surface. 
Because o[ tllis unique· 


nature of the subsurface explosions, the center of the Lagrangian arcs describing 


the a~~osphere is deliberately set at 20 m above the surface. 
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Results and Discussion 


The radius of vaporization of soil plays an important role in the 


phenomenology of the subsurface explosions. 
The structure of the blast wave 


radically changes when DOB exceeds the radius of vaporization. 
This is easily 


understood when one examines the energy coupling to the atmosphere, which is shown 


in Figure 3. 
The radius of vaporization of soil is about 2 m for 1 Kt yield, 


and the coupling efficiency rapidly decreases as DOS exceeds the radius of vapor~ 


ization. 
Figures 4, 5 and 6 clearly display the change in the structure of the 


blast waves. 
The debris boundary and the airblast for DOB = 3m appear very 


diff~rent from the others. 


In Figures 4 (DOS = 0.5 m) and 5 (OOB = 1.5 m), the shape of the Lagrangian 


grid near the ground surface suggests a Mach stem. 
There is a strong turbulence 


behind G~e wave front near the ground, and three Mach waves are observed. 
The 


wave form in Figure 6 (DOB = 3 m) is much more spherical, and only a single 


shock propagates over the ground. 


The hot source, which communicates with the ground surface, rapidly rises 


above the cold debris. 
A sphere of hot gas appears over a cone of the colder 


debris (see Figures 4 and 5). A sharp kink in the upper portion of the Lagrangian 


grid in Figures 4 and 5 is an intersection of two spherical waves. 
One wave is 


initially started by the rising ground, and the second wave is formed from the 


top of the debris when the hot gas has risen to the top. 


Figures 7 and 8 are the vertical profiles of temperature and pressure, 


respectively, along the axis of symmetry. 
They clearly show that b~e hot gas 


is expanding over the colder rarefied gas • 


The thermal radiation comes mainly from the upper portion of the debris. 


For DOS = 0.5 and 1.5 m, the source quickly penetrates the ground surface, and 


its temperature at venting is about 20 eV. The source at OOB of 3 m rises quite 
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slowly near tl1e center of the debris. 
Thus most of the radiative energy is 


absorbed by the cold materials. 
When the source finally rises above the debris, 


the temperature is down to a fraction of 1 eV. 
We examined the thermal energy 


of the debris, taking 0.25 eV to be the minimum temperature for the thermal 


spectrum. 


The thermal flux depends on the temperature and the energy density of the 


materials. 
The energy density of the radiating source (the materials at 


temperatures greater than 0.25 eV) is plotted in Figures 9, 10 and 11. 


The relative thermal yield has been measured for a number of atmospheric 
7 
explosions. 
'Ihe measured thermal yield vs the yield to mass ratio (the energy 


density) is plotted in Figure 12. 
Based on the measured thermal yield (Figure 12), 


· the computed results shown in Figures 9, 10 and 11 give the upper limits on the 


relative thermal yields of the 1 Kt source for the three DOB's. 
The upper limits 


on the thermal yield thus obtained are plotteq against tl1e depth of burial in 


Figure 13. 


'Ihere are only three points in Figure 13. 
'Ihus one can draw either a straight 


line averagin<J the three points or a dotted curve which follows the three points 


and approaches the value for a free field explosion at OOB = 0. 
The relative 


yield for the free field explosion of a 1 Kt source with the mass of 25 Kg is 


easily determined from Figure 12. It is about 15%. 


If we take the straight line in Figure 13, it implies some interesting 


results. 
The relative thermal yield of the 1 Kt source (mass 25 Kg) exploded 


at shallow depths is no more than 6%. 
The thermal yield vanishes at the depth 


of burial of about 4 m. 


When a Mach stem is formed, the shock pressure widely fluctuates near the 


ground surface. 
Above the Mach stem, the shock wave is a smoothly expanding 


spherical wave. 
Therefore, the peak overpressure of the blast wave was 


.. 


" 


5 


determined at 20 m above the ground surface and plotted against the range .in 
' 
' 


Figure 14. 
'I'he calculated overpressures are compared with the measurements 
8,9,10,11 
on Johnnie Boy, Jangle-D, and the surface bursts • 


Johnnie Boy is a 1/2 Kt device detonated at 58 em below the ground surface. 


The measured pressures are plotted against the scaled range in Figure 14. 
Jangle-D 


is a 1.2 Kt device detonated at 5.18 m below the surface. 
The computed results for 


DOB of 0.5, 1.5 and 3 m are consistent with these measurements 


A blast wave for DOB = 3 m propagates smoothly over the ground surface. 


The peak overpressures ·near the ground surface (2·m above the surface) are plotted 
12,13 
in Figure 15. 
They are compared with other calculations 
and measurements. 


1he three sets of computation done bn CEL, RAD9, and HUL~ are in good agreement 


within the computational error. 


Conclusion 


We have conducted an investigation on the effects of the subsurface nuclear 


explosions, using the current computational technology. 
For the shallow depth 


explosions, the radius of vaporization of soil plays an important role. 
The 


thermal yield and the blast strength at high pressures become a strong function 


of OOB when DOB exceeds the radius or vaporization. 


The thermal yield vanishes at DOB atout twice the radius of vaporization. 


The yield to mass ratio as well as the temperature determines the G~ermal flux. 


Since the radiative mass of the debris is almost entireiy the vaporized soil for 


n small device, the yield to mass ratio of the device itself has very little effect. 


Although DOB can change the blast strength by orders of magnitude at high 


pressures, the peak overpressures at large distances, e.g., below 10 psi, are 


affected much less by 008. The three sets of computations done on three different 


codes are mutually consistent within the computational error. 
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Closer attention should be given to the equation ·of state and .the opacity 


because the radius of vaporization is very sensitive to the equation of state, 


and the opacity is a determining factor in energy transport. 
'llle current LLL 


tables for the equation of state and the opacity are not adequate in the region 


of our interest. 
r-Dre improvement is desired in this area. 
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Figure Captions 
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Figure 3 


Figure 4 


Figure 5 


Figure 6 


Figure 7 
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